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Abstract: One of the most surprising consequences of quantum mechanics is the entan- 
glement of two or more distant particles. Although questions regarding fundamental issues 
of quantum theory still exist, quantum entanglement has started to play important roles 
in practical engineering applications. Quantum imaging is one of these exciting areas. 
Quantum imaging has demonstrated two peculiar features: (1) reproducing "ghost" im- 
ages in a "nonlocal" manner, and (2) enhancing the spatial resolution of imaging beyond 
the diffraction limit. In this article, we start with the review of classical imaging to estab- 
lish the basic concepts and formalisms of imaging. We then analyze two-photon imaging 
with particular emphasis on the physics of spatial resolution enhancement and the "ghost" 
imaging phenomenon. 



1 Introduction 

Quantum mechanics remarkably allows for the entanglement of two or more distant parti- 
cles. In a maximally entangled two-particle system, the value of an observable for each particle 
is completely undefined while the correlation between the two can be uniquely determined, 
despite the distance between them. In other words, each of the two subsystems may have com- 
pletely random values or all possible values for some physical observable in the course of their 
propagation, but the correlations of the two subsystems are determined uniquely and with cer- 
tainty through the measurement of joint-detection events p. Although questions regarding 
fundamental issues of quantum theory still exist, quantum entanglement has started to play 
important roles in practical engineering applications. Quantum imaging is one of these exciting 
areas. 

The first quantum imaging experiment was demonstrated by Pittman et al. in 1995 [2], 
inspired by the theory of Klyshko [3j . The experiment was immediately named "ghost imaging" 
due to its surprising nonlocal feature. The important physics demonstrated in that experiment, 
nevertheless, may not be the "ghost." Indeed, the original purpose of the experiment was to 
study and to test the two-particle EPR correlation in position and in momentum for an entan- 
gled two-photon system [4] [5] [6] [7] . The experiments of ghost imaging and ghost interference 
[8j together stimulated the foundation of quantum imaging in terms of geometrical and physical 
optics. 

Entangled multi-photon systems were later introduced to lithography for sub-diffraction- 
limitted imaging [9]. In 2000, Boto et. al. proposed a "noon" state and approved that the 
entangled N-photon system may improve the spatial resolution of an imaging system by a 
factor of N, despite the Rayleigh diffraction limit. The working principle of quantum lithography 
was experimentally demonstrated by DAngelo et al. in 2001 [T0| by taking advantage of an 
entangled two-photon state of spontaneous parametric down-conversion (SPDC). Due to the 
lack of a two-photon absorber, the joint-detection measurement in that experiment was on the 
Fourier transform plane rather than on the image plane. The observed Fourier transform of 
the object function is the same as the one produced by classical light of wavelength A/2. It 
was implicit in Ref.[10i that a second Fourier transform, by inserting a second lens in that 
experimental setup, would transfer the Fourier transform of the object onto its image plane, 
thus giving an image with doubled spatial resolution despite the Rayleigh diffraction limit. 
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Quantum imaging has so far demonstrated two peculiar features: (1) reproducing ghost 
images in a "nonlocal" manner, and (2) enhancing the spatial resolution of imaging beyond 
the diffraction limit. Both the nonlocal behavior observed in the ghost imaging experiment 
and the apparent violation of the uncertainty principle explored in the quantum lithography 
experiment are due to the two-photon coherent effect, which involves the superposition of 
two-photon amplitudes, a nonclassical entity corresponding to different yet indistinguishable 
alternative ways of triggering a joint-detection event, under the framework of quantum theory 
of photodetection [llj . 

The nonlocal superposition of two-photon states may never be understood classically. Clas- 
sical attempts, however, have never stopped in the history of EPR studies. Bennink et al. 
demonstrated an interesting experiment in 2002 jl2j. In that experiment, two co-rotated laser 
beams produced a projection shadow of an object mask through coincidence measurements. In- 
stead of having a superposition of a large number of two-photon probability amplitudes, Bennink 
used two correlated laser beams (imagine two back to back lasers) to simulate each two-photon 
probability one at a time. If the laser beam propagating in the object arm is blocked by the 
mask at a certain rotating angle, there would be no coincidence in that angle and consequently 
the corresponding "position" in the nonlocal "image" plane, defined by the correlated laser 
beams which were propagated to a different directions. The block-unblock of the correlated 
laser beams projected a shadow of the object mask in coincidences. Interestingly, this experi- 
ment has excited a number of discussions concerning certain historical realistic models of EPR. 
Despite the fact that the measured correlation by Bennink et al. is not the EPR correlation 
in momentum and, simultaneously, in position but rather a trivial "momentum-momentum" 
correlation defined by the two co-rotating laser beams on their Fourier transform planes, this 
experiment is a good example in distinguishing a two-photon image from a correlated projection 
shadow [6] [7]. 

Thermal light ghost imaging was another challenge. In 2004, Gatti et al. [13], Wang et 
al. [H], and Zhu et al. [T^ proposed ghost imaging by replacing entangled state with chaotic 
thermal radiation. A question about ghost imaging is then naturally raised: Is ghost imaging a 
quantum effect if it can be simulated by "classical" light? Thermal light ghost imaging is based 
on the second-order spatial correlation of thermal radiation. In fact, two-photon correlation of 
thermal radiation is not a new observation. Hanbury-Brown and Twiss (HBT) demonstrated 
the intensity spatial correlation of thermal light in 1956 16J. Differing from entangled states, 
the maximum correlation in thermal radiation is 50%, which means 33% visibility of intensity 
modulation at most. Nevertheless, thermal light is a useful candidate for ghost imaging in 
certain applications. Recently, a number of experiments successfully demonstrated certain 
interesting features of ghost imaging by using chaotic light [17] [18] [19] [20] . 

The HBT experiment was successfully interpreted as statistical correlation of intensity 
fluctuations. In HBT, the measurement is in far- field (Fourier transform plane). The measured 
two intensities have the same fluctuations while the two photodetectors receive the same mode 
and thus yield maximum correlation 

ihh) = hh + (A/1A/2). (1) 

When the two photodetector receive different modes, however, the intensities have different 
fluctuations, the measurement yields (AJ1A/2) = and gives (/1/2) = hh- One type of the 
HBT experiments explored the partial (50% ) spatial correlation (J1/2) l + 5[{pi— p2){/^9)\ of 
the thermal radiation field, where pj is the transverse coordinate of the photodetector and 

is the angular size of the source. This result has been applied in Astronomy for measuring 
the angular size of stars. 

Although Eq. ([1]) gives a reasonable explanation to the far-field HBT phenomena, the 
theory of statistical correlation of intensity fluctuation may not work for ghost imaging: (1) 
this theory fails to provide adequate interpretation for ghost imaging of entangled states. The 
visibility of ghost image of entangled two-photon state is 100% which means a 100% correlation, 
i.e., {I1I2) 5{pi ~ 92/^), where m is the magnification factor of imaging. If one insists on 
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Eq. ([T|), the mean intensities Ji and I2 must be zero. Otherwise Eq. ^ leads to non- physical 
conclusions. The measurements, however, never yield zero mean values of Ii and I2 in any 
circumstances. (2) Scarcelli et al. recently demonstrated a lens-less "near-field" ghost imaging 
of chaotic radiation [18] and pointed out that the theory of statistical correlation of intensity 
fluctuation does not work for their experiment. Differing from HBT in which the measurement 
is in far-field, Scarcelli's experiment is in near-field. In near-field, for each point on the detection 
plane, a point photodetector receives a large number of {N) modes in the measurement. The 
ratio between joint-detections triggered by "identical mode" and joint-detections triggered by 
"different modes" is N /N"^ . For a large N , the contributions from "identical mode" is negligible 
and thus (A/1A/2) = 0. Therefore, the classical idea of statistical correlation of intensity 
fluctuations will not work in the multi-mode case, as we know that different modes of chaotic 
light fluctuate randomly and independently. On the other hand, Scarcelli et al. proved a 
successful alternative interpretation based on the quantum theory of two-photon interference. 

It is interesting to see that the concept of two-photon coherence is applicable to both 
"classical" and "quantum" light. Although, two-photon superposition is a new concept bene- 
fitted from recent research of entangled states 4 , the concept is not restricted to the entangled 
states. The concept is generally true and applicable to any radiation, including "classical" ther- 
mal light. It is then reasonable to ask: What is the difference between entangled two-photon 
imaging and thermal light two-photon imaging? To what extent can two-photon correlation 
of chaotic radiation replace the roles of entangled states? Can two-photon imaging of chaotic 
light improve imaging resolution beyond the classical limit? Another reasonable question may 
trap us into an enduring debate: Is thermal light ghost imaging "classical" or "quantum"? Is 
ghost imaging, in general, "classical" or "quantum"? 

This article attempts to provide answers to these hotly debated questions. To this end, 
we start with the review of classical imaging. The treatment of classical imaging may serve to 
establish the basic concept and formalism we are going to use in the rest of this article. We 
then analyze two-photon imaging with particular emphasis on the physics of spatial resolution 
enhancement and the "ghost" imaging phenomenon. 



2 Imaging and Spatial Resolution 

The concept of imaging is well defined in classical optics. Fig. [1] schematically illustrates a 
standard imaging setup. A lens of finite size is used to image the object onto an image plane 
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Figure 1: Imaging: a lens produces an image of an object in the plane defined by the 
Gaussian thin lens equation l/s.; + I/sq ~ ^/ f- The concept of an image is based on the 
existence of a point-to-point relationship between the object plane and the image plane. 



which is defined by the "Gaussian thin lens equation" 

1 1 1 
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where So is the distance between object and lens, / is the focal length of the lens and Si is the 
distance between lens and image plane. If light always follows the laws of geometrical optics, 
the image plane and the object plane would have a perfect point-to-point correspondence, which 
means a perfect image of the object, either magnified or dcmagnified. Mathematically, a perfect 
image is the result of a convolution of the object distribution function f{po) and a (5-function. 
The (5-function characterize the perfect point-to-point relationship between the object plane 
and the image plane: 

F{Pt)^ j dpof{po)5{po + Pi/m) (3) 

J obj 

where po and pi are 2-D vectors of the transverse coordinate in the object plane and the image 
plane, respectively, and m is the magnification factor. 

Unfortunately, light behaves like a wave. The diffraction effect turns the point-to-point 
correspondence into a point-to- "spot" relationship. The (5-function in the convolution of Eq. ([3]) 
will be replaced by a point-spread function. 

F{pi)=l dpof{po)somb[ \po + Pi/m\] (4) 

J obj C 

where 

somo(x) — , 

X 

and Ji{x) is the first-order Bessel function, R is the radius of the imaging lens. The finite size 
of the spot, which is defined by the point-spread function, determines the spatial resolution 
of the imaging setup, and thus, limits the ability of making demagnified images. It is clear 
from Eq. (j4]), the use of larger imaging lens and shorter wavelength light source will result 
in a narrower point-spead function. To improve the spatial resolution, one of the efforts in 
the lithography industry is the use of shorter wavelengths. This effort is, however, limited 
to a certain level because of the lack of lenses effectively working beyond a certain "cutoff' 
wavelength. 




Figure 2: Typical imaging setup. A lens of finite size is used to produce a demagnified image 
of a object with limited spatial resolution. Replacing classical light with an entangled N- 
photon system, the spatial resolution can be improved by a factor of N, despite the Rayleigh 
diffraction limit. 

Eq. (|4]) imposes a diffraction limited spatial resolution on an imaging system while the 
aperture size of the imaging system and the wavelength of the light source are both fixed. This 
limit is fundamental in both classical optics and in quantum mechanics. Any violation would 
be considered as the violation of the uncertainty principle. 

Surprisingly, the use of quantum entangled states gives a different result: by replacing 
classical light sources in Fig. [2] with entangled N-photon states, the spatial resolution of the 
image can be improved by a factor of N, despite the Rayleigh diffraction limit. Is this a violation 
of the uncertainty principle? The answer is no! The uncertainty relation for an entangled N- 
particle system is radically different from that for N independent particles. In terms of the 
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terminology of imaging, what we have found is that the somh{x) in the convolution of Eq. ^ 
has a different form in the case of an entangled state. For example, an entangled two-photon 
system has 

X = \po + Piim\, 

So c ' 

and thus yields a twice narrower point-spread function and results in a doubling spatial resolu- 
tion for imaging. 




Figure 3: Projection: a light source illuminates an object and no image forming system is 
present, no image plane is defined, and only projections, or shadows, of the object can be 
observed. 

It should be emphasized that one must not confuse a "projection" with an image. A pro- 
jection is the shadow of an object, which is obviously different from the image of an object. 
Fig. [3] distinguishes a projection shadow from an image. The difference is more obvious when 
the projection shadow is accumulated one-ray-at-a-time, as Bennink et al. did in their exper- 
iment |12| . In a projection, the object-shadow correspondence is essentially a "momentum" 
correspondence, which is defined by the propagation direction of the light rays. This point will 
be emphasized again in the section of ghost imaging. 



3 Classical Imaging 

We now analyze classical imaging. The analysis starts from the propagation of the field 
from the object plane to the image plane. In classical optics, such propagation is described by 
an optical transfer function h{Y — ro,t — to), which accounts for the propagation of all modes 
of the field. In view of our quantum optics approach, we prefer to work with the single-mode 
propagator ^(k, r — ro,t — to) and its Fourier transform [1T][52]. It is convenient to write the 
field E{r,t) in terms of its longitudinal (z) and transverse {p) coordinates under the Fresnel 
paraxial approximation [21j : 

E{p,z,t) = J duj dn E{k,u}) g{K,uj; p, z) e~'^'^* (5) 

where E{u!, i?) is the complex amplitude of frequency u and transverse wave- vector k. In Eq. ^ 
we have taken zq = and to = at the object plane as usual. To simplify the notation, we 
have assumed one polarization. 

Based on the experimental setup of Fig. [21 g{K, oj; p, z) is found to be [2T]p2] 

g{K,uj;pi, So + Si) 

where po, pi, and pi are two-dimensional vectors defined, respectively, on the object, the lens, 
and the image planes. The first curly bracket includes the aperture function A{po) of the 
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object and the phase factor e'"'''" contributed to the object plane by each transverse mode k. 
The terms in the second and the fourth curly brackets describe free-space Fresnel propagation- 
diffraction from the source/object plane to the imaging lens, and from the imaging lens to 
the detection plane, respectively. The Fresnel propagator includes a spherical wave function 
e^f _ zk) and a Fresnel phase factor gi^lPj-PklVM^j-^k) ^ xhe third curly bracket 

adds the phase factor introduced by the imaging lens. A brief review of Fresnel diffraction- 
propagation is given in Appendix A. 

We now rewrite Eq. ([6]) into the following form 



g{K.,uj;pi, z = So + Si) 

J obi 



(27rc)^5o5i 



s 

The image plane is defined by the Gaussian thin- lens equation of Eq. Hence, the second 
integral in Eq. ([7]) simplifies and gives, for a finite sized lens of radius R, the so called point- 
spread function of the imaging system: somb{x) = 2Ji{x)/x, where x = [j- ^\po + Pi/'m\\, Ji[x) 
is the first-order Bessel function and m ~ Si/so is the magnification of the imaging system. 

Substituting the result of Eqs. ([7]) into Eq. ([5]) enables one to obtain the classical self- 
correlation of the field, or, equivalently, the intensity on the image plane 

I{pi,Zi,U) = {E* {p^, z^,U) E{pi, Zi,ti)) (8) 

where (...) denotes an ensemble average. We assume monochromatic light for classical imaging 
as usual ^1] . 

Case (I): incoherent imaging. The ensemble average of {E*{k,lu) E{k' ,uj)) yields zeros 
except when k, ~ k' . The image is thus 

lip,) oc f dp,\A{p,)\^\somb[--\po + ^\]\'. (9) 
J ' ' ' So c m ' 

An incoherent image, magnified by a factor of to, is thus given by the convolution between 
the squared moduli of the object aperture function and the point-spread function. The spatial 
resolution of the image is thus determined by the finite width of the |soTO5p-function. 

Case (II): coherent imaging. The coherent superposition of the k, modes in both E*{pi,T) 
and E{pi,T) results in a wavepacket. The image, or the intensity distribution on the image 
plane, is thus 

/(ft)oc / dpoAipo)e'^\P"\%omb[--\po + -\]^- (10) 

Jobj So C TO 

A coherent image, magnified by a factor of to, is thus given by the squared modulus of the 
convolution between the object aperture function (multiplied by a Fresnel phase factor) and 
the point-spread function. 

For Si < So and So > /, both Eqs. ^ and PH)) describe a real demagnified inverted image. 
In both cases, a narrower soTO6-function yields a higher spatial resolution. Thus, the use of 
shorter wavelengths allows for the improving of the spatial resolution of an imaging system. 



4 Two-photon Imaging 

To demonstrate the working principle of quantum imaging, we replace the classical light 
source in Fig. [3 with an entangled two-photon source such as spontaneous parametric down- 
conversion (SPDC) [53] [3] and replace the ordinary film with a two-photon absorber, which is 
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sensitive to two-photon transition only, on the image plane. We will show that, in the same 
experimental setup of Fig. [21 an entangled two-photon system gives rise, on a two-photon 
absorber, to a point-spread function twice narrower than the one obtained in classical imaging 
at the same wavelength. Then, without employing shorter wavelengths, entangled two-photon 
states improve the spatial resolution of a two-photon image by a factor of 2. 

What is special about the entangled two-photon state of SPDC? Let us imagine a mea- 
surement in which we place two point-like single photon counting detectors (Z?i and D2) on 
the output surface of an SPDC crystal. The nearly coUinear signal-idler system generated by 
SPDC can be described by the two-photon state: 

*) = *o dKsdKt5{Ks + Ki) I (hjsdi^i5{uJs+Ui-Up)a\ns,i^s)a\K.i,uji)\Q) (11) 



where tOj , (j = s,i,p), is the frequency and the transverse wavevector of the signal, idler, and 
pump. To simplify the calculation, we have assumed a CW single frequency plane-wave pump 
with Kp = 0. The joint-detection probability between the point-like detectors Di and D2 located 
at (yOi,zi) and (/92,-Z2), respectively, is calculated from the Glauber theory of photodetection 

G(2)(p'i,Zi,ti;p2,Z2,t2) = \{0\Ei+\p2,Z2,t2)E[+\pi,Z,,t,)\M')\^ 

= \-^{pi,ZiM]P2,Z2M)?- (12) 

At this point we are taking zi = Z2 = 0. The transverse part of the effective two-photon 
wavefunction ^'(pi,p2) is thus calculated as: 

■^{pi,P2)~-^o5{pi-p2)- (13) 

Both Eqs. (fTT|) and suggest that the entangled signal-idler photon pair is characterized 
by EPR correlation [Tj in transverse momentum and transverse position; hence, similar to the 
original EPR state, we have [7]: 

A(k, + = & A(pi - P2) = (14) 
with Aks ~ cx), A/ii ~ 00, Api ^ 00, A/92 ^ 00. 

In EPR's language, the signal and the idler may come out from any point on the output plane 
of the SPDC. However, if the signal (idler) is found in a certain position, the idler (signal) must 
be found in the same position, with certainty (100%). Furthermore, the signal and the idler 
may have any transverse momentum. However, if a certain value and direction of the transverse 
momentum of the signal (idler) is measured, the transverse momentum of the idler (signal) will 
certainly have equal value and opposite direction. 

Concerning the experimental setup of Fig. ^ in light of Eq. (1141) . if the object is placed 
very close to the output plane of the SPDC crystal, the signal photon and the idler photon are 
guaranteed to come out from the same point of the object plane {ps = Pi) and to stop at the 
corresponding unique point on the image plane, thus forming an image together (i.e., a two- 
photon image), whether enlarged or demagnified. In addition, the signal- idler pair propagates 
in a special way in which the transverse momenta of the pair must satisfy Ks + K,i = and 
thus undergo "two-photon diffraction." Two-photon diffraction is radically different from that 
of two independent photons and yields a different spatial resolution of the image. This unique 
feature will be shown in the following calculation. 

The fields at the detectors are the same as in the classical imaging case. By inserting the field 
operators into \E'(/9i, zi, ii; /92, ^2, ^2), as defined in Eq. (|12[) . and considering the commutation 
relations for the creation and the annihilation operators, we find for the effective two-photon 
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wavefunction in the image plane 



'^{Pl,Zl,ti;P2,Z2,t2) (15) 

. *o / d,. dU. sr.. + n.) J d.. d.. + - 

X g{Ks,u!s;pi,zi) e~*"=*i g(Kj, uj.^; p2, 22) e"*"**^ 

Substituting the g functions of Eq. ([7]) into Eq. (jlSp and applying the (5-functions in the 
state to reduce the variables of the integral in Eq. (fT5|) from 4 to 2, we obtain 

*(pi,n;p2,T-2) 

(X y dpo A(po) e^^l'^°l' j dp', A{p',) e'^l'^-'ol' 

X e"'5^l'^°l'e*5^'''''°l'{y"dAt,e'^'-['^°-''"*'°l} (16) 

where Tj = tj — Zj/ c, v is defined from = + v and w^. = i-Lip/2 — following a;^ + cjj = cjp. 

Now we consider a two-photon film (two-photon absorber) , or equivalently scan Z?i and D2 
together, on the image plane to achieve pi = P2, zi — Z2 and examine the two integrals in the 
two curly brackets. It is easy to see that the integral of dK. yields S{po — p' o), which is consistent 
with Eq. (|13p . The integral of di' gives a similar (5-function in the form of S[{-^ + -^){pi — p' i)] 

while taking po — p' Pi = P2 = P, and ri = T2. These results indicate that the propagation- 
diffraction of the signal photon and the idler photon are not independent. The "two-photon 
diffraction" couples the two integrals in po and p', as well as the two integrals in pi and p'l and 
gives the G'^^-' function 



y So c 



G(2)(p,p)oc / dpoA^(p-)e-^l^^°l ^ ' ^ (17) 



obj 



\s„ c \P° ^ m 



which indicates that a coherent image magnified by a factor oi m ~ Si/ Sq is reproduced on the 
image plane by joint-detection or by two-photon absorption. 

In Eq. (jl7p . the point-spread function is characterized by the pump wavelength Ap = As^i/2; 
hence, the point-spread function is twice narrower than in the (first order) classical case. An 
entangled two-photon state thus gives an image in joint-detection with double spatial resolution 
when compared to what one would obtain in classical imaging. Moreover, the spatial resolution 
of the two-photon image obtained by SPDC is further improved because it is determined by 
the function 2Ji{x)/x^, which is much narrower than the somb{x). 

It is interesting to see that, different from the classical case, the frequency integral over 
Aujs 00 does not give any blurring problem, but rather enhances the spatial resolution of 
the two-photon image. These characteristics of entangled states may be useful for certain 
applications ^4J. 

Can we replace entangled states with classical light in the same setup of Fig. [2] and still 
have a twice narrower point-spread function in two-photon joint-detection? 

We will first quickly examine the case of coherent radiation. By definition, coherent light 
is characterized by the relation: G'-^-' = jG'-^-'p [11;. The two-photon image produced by laser 
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light is readily obtained from Eq. ([TO)) : 



dpo A{po 



obj 



dpo A{po) e 



somO \Po^ — 

So c m 



obj 



somby |po H 1 

So c m 



(18) 



which is simply a product of two independent first-order coherent images (see Eq. (|10l) ). 

Now we analyze chaotic light. The standard form of the G^^^-function of chaotic light is 



G'^'\p,,P2) - g(V(/?i)gW(p2) + \g[\\p,,P2)\ 



(19) 



where G^\^(pi) and G^22 iP'^) mean intensity distributions on the image plane. It is 

straightforward to find that the first term of Eq. p^. i.e., gJ\' (pi) G22^ iP2), is a simple product 
of two independent incoherent classical (first-order) images. The interesting part of Eq. ([TO]) is 
in the second term where 



Gi2^(pi,2i,ii;P2,2:2,i2) / dfig*{K,uj;pi,zi)e'"^*^ g{K,uj;p2,Z2)e~ 



(20) 



Substituting the g* and g functions into Eq. (j20p and considering Di and D2 on the image 
plane, we obtain 



g[^^{Pi,P2) oc 



obj 



dpoA*{po) e *2c.. 



Ip'o 



obj 



dp', e 



(21) 



Similar to the first-order imaging, in Eqs. ([20]) and ([21]) we have assumed monochromatic light 
to avoid "aberration" complications. Differing from that of SPDC, here an integral of Auj ^ 00 
results in a constant G^2'' , which is consistent with our experimental observation (17j . 
Evaluating the integrals in Eq. ([2T|) . the G^^GjV term turns out to be 



\g[]^{Pi,P2)\'^ oc I dpo\A{po)f somb(—- 

J obj ^So C 



Po H 

m 



\ somb( 



Po H 

m 



(22) 



Unlike the G^^^Ga^ term, here, we do not have two independent images. This is the reason 
that chaotic light has been successfully used to mimic certain nonlocal features of entangled 
systems [18]. Chaotic light, however, is not a viable alternative to entangled sources for quantum 
lithography. In fact, when Di and D2 are scanned together, or a two-photon sensitive material 
is employed for recording the photon pair, i.e., pi = P2 = p, Eg. 1221 becomes: 



\G['Jip,p)\^ (X I / dpo \Aipo)\^ \somb[ 

J obj 



RUl: 



Po 



(23) 



which is the same result as given by the G^^-'Gj^ term. Namely, it is simply the product of 
two independent incoherent first-order classical images. This conclusion relates quite nicely to 
the findings of 25J: in that experiment, a doubly modulated interference-diffraction pattern 
was observed by scanning the detectors in opposite directions on the Fourier transform plane. 
However, no second-order interference pattern can be observed if the detectors are scanned 
together. 
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5 Ghost Imaging 



Pittman et al demonstrated a surprising two-photon imaging experiment in 1995 The 
experiment was immediately named "ghost imaging" due to its "nonlocal" feature. The impor- 
tant physics demonstrated in the experiment, nevertheless, may not be the so called "ghost." 
Indeed, the original purpose of the experiment was to study the EPR correlation in position 
and in momentum and to test the EPR inequality [H [71 [S] for the entangled signal-idler photon 
pair of SPDC. 
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Figure 4: Schematic set-up of the "ghost" image experiment. 

The schematic setup of the "ghost" imaging experiment is shown in Fig. [H A CW laser 
is used to pump a nonlinear crystal, which is cut for degenerate type-II phase matching to 
produce a pair of orthogonally polarized signal (e-ray of the crystal) and idler (o-ray of the 
crystal) photon. The pair emerges from the crystal coUinear, with uig LOi = The pump 

is then separated from the signal-idler pair by a dispersion prism, and the remaining signal and 
idler beams are sent in different directions by a polarization beam splitting Thompson prism. 
The signal beam passes through a convex lens with a AQQmm focal length and illuminates a 
chosen aperture (mask). As an example, one of the demonstrations used letters "UMBC" for 
the object mask. Behind the aperture is the "bucket" detector package Di, which consists 
of a short focal length collection lens in whose focal spot is an avalanche photodiode. Di is 
mounted in a fixed position during the experiment. The idler beam is met by detector package 
I?2, which consists of an optical fiber whose output is mated with another avalanche photodiode. 
The input tip of the fiber is scanned in the transverse plane by two step motors. The output 
pulses of each detector, which are operating in photon counting mode, are sent to a coincidence 
counting circuit for the signal-idler joint-detection. 

By recording the coincidence counts as a function of the fiber tip's transverse plane coor- 
dinates, the image of the chosen aperture (for example, "UMBC") is observed, as reported in 
Fig. O It is interesting to note that while the size of the "UMBC" aperture inserted in the 
signal beam is only about 3.5mm x 7mm, the observed image measures 7mm x 14mm. The 
image is therefore magnified by a factor of 2. The observation also confirms that the focal 
length of the imaging lens /, the aperture's optical distance from the lens So, and the image's 
optical distance from the lens Si (which is from the imaging lens going backward along the 
signal photon path to the two-photon source of SPDC crystal then going forward along the 
path of idler photon to the image), satisfy the Gaussian thin lens equation. In this experiment, 
So was chosen to be So — 600mm, and the twice magnified clear image was found when the 
fiber tip was on the plane of Si = 1200mm. When D2 was scanned on transverse planes not 
determined by the Gaussian thin lens equation the images blurred out. 
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Figure 5: (a) A reproduction of the actual aperature "UMBC" placed in the signal beam, 
(b) The image of "UMBC" : coincidence counts as a function of the fiber tip's transverse 
coordinates in the image plane. The step size is 0.25mm. The image shown is a "slice" at 
the half maximum value. 

The measurement of the signal and the idler subsystem themselves are very different. The 
single photon counting rate of £>2 was recorded during the scanning of the image and was found 
fairly constant in the entire region of the image. This means that the transverse coordinate 
uncertainty of either signal or idler is considerably large compared to that of the transverse 
correlation of the entangled signal-idler photon pair: Axi (Ayi) and Ax2 (Ay2) are much 
greater than A{xi — X2) (A(yi — 2/2))- 




S„ Si 



Figure 6: An unfolded setup of the ghost imaging experiment, which is helpful for under- 
standing the physics. Since the two-photon "light" propagates along "straight-lines", it is 
not difficult to find that any geometrical light point on the object plane corresponds to an 
unique geometrical light point on the image plane. Thus, a ghost image of the object is 
made nonlocally in the image plane. Although the placement of the lens, the object, and 
detector D2 obeys the Gaussian thin lens equation, it is important to remember that the 
geometric rays in the figure actually represent the two-photon amplitudes of an entangled 
photon pair. The point to point correspondence is the result of the superposition of these 
two-photon amplitudes. 

The EPR (5-functions, S{ps — Pi) and 6{ks + Ki) in transverse dimension, are the key to 
understand this interesting phenomenon. In degenerate SPDC, although the signal-idler photon 
pair has equal probability to be emitted from any points on the output surface of the nonlinear 
crystal, the transverse position ^-function indicates that if one of them is observed at one 
position, the other one must be found at the same position. In other words, the pair is always 
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emitted from the same point on the output plane of the two-photon source. Simultaneously, 
the transverse momentum (5-function defines the angular correlation of the signal-idler pair. 
The transverse momenta of a signal-idler amplitude are always equal but pointed in opposite 
directions, which means that the two-photon amplitudes are always existing at 

roughly equal yet opposite angles relative to the pump. This then allows for a simple explanation 
of the experiment in terms of "usual" geometrical optics in the following manner: we envision 
the nonlinear crystal as a "hinge point" and "unfold" the schematic of Fig. 2] into the Klyshko 
picture of Fig. [6l The signal-idler two-photon amplitudes can then be represented by straight 
lines (but keep in mind the different propagation directions) and therefore, the image is well 
produced in coincidences when the aperture, lens, and fiber tip are located according to the 
Gaussian thin lens equation of Eq. Q . The image is exactly the same as one would observe on 
a screen placed at the fiber tip if detector Di were replaced by a point-like light source and the 
nonlinear crystal by a reflecting mirror. 

Following a similar analysis in geometric optics, it is not difficult to find that any geometrical 
"light spot" on the object plane, which is the intersection point of all possible two-photon 
amplitudes coming from the two-photon light source, corresponds to an unique geometrical 
"light spot" on the image plane, which is another intersection point of all the possible two- 
photon amplitudes. This point to point correspondence made the "ghost" image of the subject- 
aperture possible. Despite the completely different physics from classical geometrical optics, the 
remarkable feature is that the relationship between the focal length of the lens /, the aperture's 
optical distance from the lens So^ and the image's optical distance from the lens S'i, satisfy 
the Gaussian thin lens equation of Eq. Although the placement of the lens, the object, 
and the detector D2 obeys the Gaussian thin lens equation, it is important to remember that 
the geometric rays in the figure actually represent the two-photon amplitudes of a signal-idler 
photon pair and the point to point correspondence is the result of the superposition of these two- 
photon amplitudes. The "ghost" image is a realization of the 1935 EPR gedankenexperiment. 

Now we calculate G^'^\po,pi) for the "ghost" imaging experiment, where po and pi are the 
transverse coordinates on the object plane and the image plane. We will show that there exists 
a 5-function like point-to-point correlation between the object plane and the image plane, i.e., 
if one measures the signal photon at a position of po on the object plane the idler photon can 
be found only at a certain unique position of pi on the image plane satisfying 5{po — pi/m), 
where m = —{si/so) is the image-object magnification factor. After approving the (5-function 
correlation, we show how the object function of A{po) is transferred to the image plane as a 
magnified image A{pi/m). Before starting the calculation, it is worth to emphasize again that 
the "straight lines" in Fig. [6] schematically represent the two-photon amplitudes all belong to a 
pair of signal-idler photon. A "click-click" joint measurement at (ri, ii), which is on the object 
plane, and (r2,t2), which is on the image plane, in the form of EPR ^-function, is the result of 
the coherent superposition of all these two-photon amplitudes. 

We follow the unfolded experimental setup shown in Fig.[7]to establish the Green's functions 
g{Ks,^s, Po, Zo) and ^(Ki, w.;, p2, 22)- In arm-1, the signal propagates freely over a distance di 
from the output plane of the source to the imaging lens, then passes an object aperture at 
distance So, and then is focused onto photon counting detector Di by a collection lens. We 
will evaluate g{Ks, ^s, Po, Zq) by propagating the field from the output plane of the two-photon 
source to the object plane. In arm-2, the idler propagates freely over a distance d2 from the 
output plane of the two-photon source to a point-like detector D2. g{^i, oJi, P2, 22) is thus a free 
propagator. 

(I) Arni-1 (source to object): 

The optical transfer function or Green's function in arm-1, which propagates the field from 
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Figure 7: In arm-1, the signal propagates freely over a distance di from the output plane of 
the source to the imaging lens, then passes an object aperture at distance So, and then is 
focused onto photon counting detector Di by a collection lens. In arm-2, the idler propagates 
freely over a distance d2 from the output plane of the source to a point-like photon counting 
detector Z?2. 
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the source plane to the object plane, is given by: 
g{R.s,uJs;Po, Zo^di + So) 
= e'^''" / dpi / dps { 

J lens J source 

X e-^wl'^'l' l^^e^^!'''-''-!'), (24) 
L 27rcSo J 

where ps and pi are the transverse vectors defined, respectively, on the output plane of the source 
and on the plane of the imaging lens. The terms in the first and second curly brackets in Eq. (|24p 
describe free space propagation from the output plane of the source to the imaging lens and 
from the imaging lens to the object plane, respectively. Again, e^^^l p=-P' I ^nd e'^cso I 
are the Fresnel phases we have defined in Appendix A. Here, we treat the imaging lens as a 
thin- lens. Thus, the transformation function of the imaging lens is approximated as a Gaussian 



(II) Arm-2 (from source to image): 

In arm-2, the idler propagates freely from the source to the plane of D2, which is also the 
plane of the image. The Green's function is thus: 



g{^,,ujf,p2,Z2 = d2)^—-^e'^''' dp'^e'^^"^-"'^ e'^'-c^ (25) 

27rcd2 Jsource 

where p'^ and p2 are the transverse vectors defined, respectively, on the output plane of the 
source, and on the plane of the photodetector Z?2- 

(III) ^{po,pi) and G^'^\poT Pi) (object plane - image plane): 

To simplify the calculation and to focus on the transverse correlation, in the following 
calculation we assume degenerate {lOs = LOi — oj) and coUinear SPDC. The transverse two- 
photon effective wavefunction '^{po, P2) is then evaluated by substituting the Green's functions 
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g{K,s, po, Zo) and g{K,i, uj; p2, ^2) into the expression given in Eq. (|15p : 

(X e'^^^""'^'''^ J diis dKi6{Ks + Ki) 

dpi [ dp; e*'""-'^=e'5S3rl'^=>-w I' 



X 



ens source 



X e 



IPil^ e*2^'^'~^° 



X / dp, e'^^-P'^ e"^^"''-'^'^ (26) 

source 

where we have ignored all the proportional constants. Completing the double integral of dn, 
and dKs 

J dKs dK, SiKs + e"'^'-P' e'^^-P'^ ^ 5[ps - p1), (27) 

Eq. (|26|) becomes: 

*(Po,P2) (28) 

lens J source 

We then complete the integral on dps, 

^'(Po,/02) oce*'^("«+'^') I dpje*5S-l'^^"'^'l'e~'wlP'l'e'2^l^-'-^-°l', (29) 

J lens 

where we have replaced di + d2 with Si (as depicted in Fig. [7|). Although the signal and 
idler propagate to different directions along two optical arms, Interestingly, the Green function 
in Eq. (j29p is equivalent to that of a classical imaging setup, if we imagine the fields start 
propagating from a point po on the object plane to the lens and then stop at point p2 on the 
imaging plane {p2 ^ pi)- The mathematics is consistent with our previous qualitative analysis 
of the experiment. 

The integral on dpi yields a point-to-point relationship between the object plane and the 
image plane that is defined by the Gaussian thin-lens equation: 

dpie'*[^+*-7llp-|%-'^T(lt+lt)-p- + (30) 



where we have replaced p2 with pi. In Eq. (|30p . the integral is approximated to infinity and 
the Gaussian thin- lens equation of Eq ([2]) is apphed. We have also defined m — Si/so as the 
magnification factor of the imaging system. The function 5{po + Pijrri) indicates that a point 
of Po on the object plane corresponds to a unique point of pi on the image plane. The two 
vectors pointed to opposite directions and the magnitudes of the two vectors hold a ratio of 

= \Pi\l\Po\- 

If the finite size of the imaging lens has to be taken into account (finite radius i?), the 
integral yields a point-spread function of somb{x): 

[ dpie^'^^^ + ^^-P' o,somb(^'^[po + ^]) (31) 

J lens ^So C m / 

where, again, somb{x) — 2Ji{x)/x, Ji{x) is the first-order Bessel function. The point-spread 
function turns the point-to-point correspondence between the object plane and the image plane 
into a point-to- "spot" relationship and thus limits the spatial resolution. 
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Therefore, by imposing the condition of the Gaussian thin-lens equation, the transverse 
two-photon effective wavefunction is approximated as a (5 function 

*(Po, ft) - 5{po + Pi/m), (32) 

which indicates a point to point EPR correlation between the object plane and the image plane, 
i.e., if one observes the signal photon at a position of po on the object plane, the idler photon 
can only be found at a certain unique position of pi on the image plane satisfying po + Pi/m = 
with m — Si/ So- 

We now include an object-aperture function, a collection lens and a photon counting de- 
tector Di into the optical transfer function of arm-1 as shown in Fig.[4l The collection-lens- I?i 
package can be simply treated as a "bucket" detector. The "bucket" detector integrates all 
^'(/9o,P2) that pass the object aperture A{po) as a joint photodetection event. This process is 
equivalent to the following convolution: 

i?i,2 0c/ dpo\A{po)\^\^l{po,p^)\^ ■^\A{p,/m)\^ (33) 

object 

where, again, D2 is scanning in the image plane, p2 = ft- 

As we have discussed earlier, the position-position EPR correlation is the result of the 
coherent superposition of two-photon amplitudes. 




In principle, one signal-idler pair contains all the necessary two-photon amplitudes that generate 
the point-to-point correspondence between the object and the ghost image plane. To emphasize 
this concept, we name this kind of ghost image as two-photon coherent image. 

We emphasize again, one should not confuse a classical "momentum-momentum" correla- 
tion with the EPR correlation in position and in momentum. Figure [5] is a schematic picture 
of the experiment of Bennink et al. [12j . which distinguishes a trivial classical momentum- 
momentum correlation from EPR. 




Lens-1 Lens-2 



Figure 8: The projection shadow of an object mask is reproduced in coincidences between 
Di and D2 by the use of two co-rotating correlated laser beams. The object-shadow cor- 
respondence is an accumulation of the "momentum-momentum" correlation between each 
pair of "shots" . 



6 Ghost Imaging of Chaotic Light 

Unlike first-order correlation, which is considered as a coherent effect of the electromagnetic 
field, the second-order correlation of radiation is usually considered as the statistical correlation 
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of intensity fluctuations. The first set of second-order correlation of thermal light was demon- 
strated in 1956 by Hanbury Brown and Twiss (HBT) with two different type of correlations: 
temporal and spatial [16] . The HBT experiment created quite a surprise in the physics commu- 
nity with an enduring debate about the classical or quantum nature of the phenomenon. It has 
been popular to consider that the HBT experiment measures the classical statistical correlation 
of the intensity fluctuations of the radiation [26] : 



( A/i A/2 ) = ((/i - h){h - 12)) = ( hh ) - /1/2 



(35) 



where Ii and I2 are the mean intensities of the radiation measured by photodetectors Di and 
D2, respectively. 



(a) HBT -Far Field 



(b) Ghost Image - Near Field 



I 




Figure 9: (a) Hanbury Brown and Twiss configuration, (b) lens-less Ghost imaging config- 
uration. 



Figure [9] (a) is a schematic of the historical HBT experiment which measures the second- 
order transverse spatial correlation of a monochromatic radiation of wavelength A coming from 
a distant star with an angular size of /S.9. Both photodetectors Di and D2 can be scanned along 
the axes xi and X2^ here, we assume 1-D scanning to simplify the discussion. The second-order 
transverse spatial correlation function r''^)(a;i, X2) is expected to be 

' {xi,X2) = {hh) ^ + smc [ J | (36) 

where we have simplified the problem to 1-D and assumed h = h = h- The second term in 
Eq. ([M]) . /q sin(?['K/S.9{xi + X2)/A], is interpreted as the correlation of intensity fluctuations. 
This term is useful in astronomy for angular size measurement of stars. For short wavelengths, 
this function quickly drops from its maximum to minimum when xi + X2 goes from zero to a 
value such that l^9[xi -f X2)/A = 1. Thus, we effectively have a "point" to "point" relationship 
between the xi and X2 plane, i.e., for each positive (negative) value of xi there exist only one 
negative (positive) value of X2 that may have nonzero correlation of intensity fluctuation. In 
fact, the planes of xi and X2 are in the far-field zone of the finite-size distant star (equivalent 
to the Fouier transform plane). Therefore, the measured quantity is the correlation between 
the transverse k vectors of the radiation. For a narrow function, the non-zero correlation 
corresponds to the case of equal transverse wavevectors: ki = K2- This is consistent with the 
physics behind the model of classical correlation. It is natural to imagine that the radiation 
coming from the same mode of the electromagnetic field, passing through the same optical 



16 



path, would have identical intensity fluctuations, while radiation coming from different modes, 
passing through different optical paths would not share the same intensity ffuctuations. 

Quantum models of HBT experiment derives the same correlation function [27] [28] . How- 
ever, if "classical" works, why quantum? The classical statistical interpretation has been widely 
accepted. Moreover, the concept of intensity fluctuation correlation has been even extended to 
quantum models to take over the concept of two-photon coherence. The philosophy of "pho- 
ton bunching" is essentially a phenomenological extension to quantum theory of the statistical 
correlation on photon number fluctuations. 

In the past twenty years, the massive research on quantum entanglement has brought new 
challenges to the classical statistical correlation interpretation. For example, replacing the 
chaotic light with an EPR type two-photon entangled state in Fig. [§] (a), the second-order 
transverse spatial correlation function turns out to be 

{hh} ^ Iq sine [ J. (37) 

Based on the concept of classical statistical correlation of intensity fluctuation, the mean inten- 
sities Ii and I2 must be zero in this case, otherwise Eq. ([35]) leads to non-physical conclusions. 
The measurements, however, never yield zero mean values of Ii and I2 in any circumstances. 

Thus the concept of classical statistical correlation of intensity fluctuation may not work for 
entangled two-photon states. Two-photon correlation experiments with entangled photons have 
been explained in terms of the superposition of indistinguishable alternatives, i.e., two-photon 
probability amplitudes, that can lead to a joint-detection event [4]. Such alternatives, however, 
represent a troubling concept in classical theories, because this concept has no counterpart in 
classical electromagnetic theory of light and it is nonlocal. If accepted, the nonlocal behavior 
of the radiation has been classified as a peculiar property of non-classical sources. 

More interestingly, we further ask ourselves: does the statistical correlation of intensity 
fluctuation always work for chaotic light? We are going to analyze a recent ghost imaging 
experiment of chaotic light by Scarcelli et ai, aimed at answering this question: "Can two- 
photon correlation of chaotic light be considered as correlation of intensity fluctuations?" [TH] 
We will conclude that the second-order correlation of chaotic light is a two-photon interference 
phenomenon. 

Figure[n](b) illustrates the setup of the experiment. Radiation from a chaotic pseudothermal 
source |29| was divided in two optical paths by a non-polarizing beam splitter. In arm A an 
object, a double slit with slit separation d = 1.5mm and slit width a = 0.2mm, was placed at 
a distance dA = 139mm and a bucket detector (Di) was just behind the object. In arm B a 
point detector D2 scanned the transverse planes in the neighborhood of a distance ds = dA 
from the source. The correlation was measured by either photon-counting-coincidence circuit 
or by standard HBT correlator. In the photon counting regime, two Geiger mode avalanche 
photodiodes were employed for single-photon level measurement. In HBT scheme, the two 
photodetectors were Silicon PIN photodiodes. The bucket detector Di was simulated by using 
a short focal length lens (/ = 25mm) to focus the light onto the active area of the detector 
while the point detector D2 was simulated by a pinhole. 

Figure [TU] reports the measured two-photon image of the double-slit. The result shows an 
equal-size reproduction of the double slit when scanning photodetector D2 along X2 axis, which 
is located at distance ds — dA from the source. In Fig. [TUl the DC background is subtracted 
from the image. Here, again, we have made a 1-D scanning to simplify the discussion. 

Let us flrst clarify the main differences of this measurement compared to Hanbury Brown 
and Twiss types of experiments. As we described in the beginning of this section, the measure- 
ments of HBT are in the far-field zone, which measure the momentum-momentum correlation 
of the field. In the reported experiment, instead, the measurements were in the near field zone 
{A0 ^ lOX/d) and therefore effectively measured the position-position correlation between the 
object plane and the image plane. The observation is a lensless two-photon ghost image. 
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Figure 10: The results of Lensless imaging. Normalized two-photon correlation vs transverse 
position of D2. The DC constant is subtracted from the correlation. 




Figure 11: Schematic experimental setup of a secondary image measurement. 



To confirm that the observed results corresponded to an image, Scarcelli et al. constructed 
a secondary imaging system, which is illustrated in Fig. [TTl by using a f — 85mm lens to 
image the ghost image onto a secondary image plane. The imaging lens was located at a 
distance of ds ~ 253mm from the source. A magnified secondary image was observed at 
d'g ^ 330mm from the lens, which is given by the Gaussian thin lens equation, by scanning 
D2 on the transverse plane. Fig. [12] reports the measured magnified secondary image of the 
ghost image with the expected magnification M = d'g/{dB — dA) ^ 2.9. For this examination 
this experiment used two masks with more complicated structures: one. Fig. I12r a). with the 
starting letter of the cities of the authors (2.3 x 2.5 mm) and the other. Fig. [T2r b). with the 
acronym of their institution (6 x 1.9 mm). Fig. ll2r c) shows the image obtained with the actual 
revealed correlation measurements to show the high contrast of the image, while in Fig. [T2yd) 
is shown a "slice" of the image around the half maximum of the correlation. 

The lens-less ghost imaging setup of Scarcelli et al. is a straightforward modification of 
the HBT. One needs simply move the two HBT photodetector from far-field to near-field. We 
cannot stop to ask: What has been preventing this simple move for 50 years (1956-2006)? 
Something must be terribly misleading to give us such misled confidence not even to try the 
near-field measurement in half a century. 

The classical model in terms of statistical correlation of intensity fluctuations would not 
work for Scarcelli's experiment. Unlike the HBT experiment in which the correlation was mea- 
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I 17.2 mm | 

Figure 12: Results of the secondary Imaging. On the left, (a) and (b) are the two object 
used, on the right are shown the obtained images: (c) is a density plot of the normalized 
correlation, (d) is a "slice" of the measured correlation around the half maximum. 

sured in the far-field zone. Scarcelli's experiment is in near-field. In the near-field configuration, 
any point on the object plane and on the image plane are "hit" by many different k vectors 
(modes). Fig. [T3l schematically illustrates the situation. Assuming a chaotic radiation source of 
finite dimension consists of a large number of N point sub-sources. By definition of "chaotic", 
each point sub-sources radiate and fluctuate independently and randomly. Now, we place two 
point photodetectors in near-field as illustrated in Fig. [13] and make a similar correlation mea- 
surement as that of HBT. It is easy to see that the chance of receiving radiation from the same 
sub-source is less then that of receiving from different sub-sources. The ratio between the two 
cases is approximately N/N'^. For a large number of N, the contributions from "identical sub- 
source" is negligible and thus (A/1A/2) — 0, as we know that different sub-sources of chaotic 
light fluctuate randomly and independently. HBT is a special case. The two photodetectors 
are placed in far-flcld. A far-field plane is equivalent to the Fourier transform plane. On the 
Fourier transform plane, each point corresponds to a mode of radiation, and corresponds to 
a point sub-source. Thus, the chance for the two photodetectors to receive "identical mode" , 
or receive radiations from "identical point sub-sources" , is non-negligible. In this situation, 
(A/1A/2) ^ 0. 




Figure 13: A simplified model of near-field joint detection of chaotic source. It is easy to see 
that the chance of receiving radiation from the same sub-source is less then that of receiving 
from different sub-sources. 

On the other hand, the quantum model for Scarcelli's experiment is straightforward. Even 
if the source is a "classical" chaotic light, the ghost image is explained as a quantum two-photon 
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effect. In a joint photodetection event of chaotic radiation, due to the superposition of two- 
photon ampHtudes, there exists 50% chance for a precise position-position correlation, i.e., if 
photon 1 is measured at pi, photon two can only be measured at a unique certain position of 
P2 ■ The following is a brief quantum mechanical calculation for the second-order correlation of 
chaotic radiation. 

In the quantum theory of photodetection , the second-order correlation function is given 
in Eq. What we need to do now is to find the state of the radiation and the operators of 

the field. In the photon counting regime, it is reasonable to model the thermal light in terms of 
single photon states. We have provided a simple model of thermal state in Appendix B. Similar 
to our earlier discussion, we will focus on the transverse correlation and assume the thermal 
radiation monochromatic: 

p ^\Q){Q\ + \et^^a\K)a) {K')\{)){Q\a{ti')a{K) (38) 

where |e| 1. Basically we are modeling the light source as an incoherent statistical mixture 
of two photons with equal probability of having any transverse momentum k and k! . 
The transverse spatial part of the second order correlation function is thus: 

= Y^\{Q\e'^+\p,)e[+\p,)\UI^^,)\\ (39) 

where pj is the transverse coordinate of the j*'* detector. The transverse part of the electric 
field operator can be written as follows: 

a 

where a(K) is the annihilation operator for the mode corresponding to k and gjiPj',K) is the 
Green's function associated to the propagation of the field from the source to the j*'' detector 

m- 

Substituting the field operators into Eq. we obtain: 

G^^Hpi;P2) = E \g2iP2,K)gi{pi,K') + g2{P2,K')gi{Pl,K)\'^- (41) 
k,k' 

This expression represents the key result toward the understanding of the phenomenon. In fact, 
it expresses an interference between two alternatives, different yet equivalent, which lead to a 
joint photodetection. The interference phenomenon is not, as in classical optics, due to the su- 
perposition of electromagnetic fields at a local point of space-time. It is due to the superposition 
of 52(^2, K)gi{pi, k') and (72 (p2, K')giipi, k), the so-called two-photon amplitudes, non-classical 
entities that involve both arms of the optical setup as well as two distant photodetection events. 
Eq. ([4T|) can be further evaluated into the form of 

G(2)(Pi;P2) oc E|.9l(pi,K)|2E|g2(/?2,K')P + |^5l*(Pl,^)52(P2,K)|^ 

^ G[\\p,)Gi\\p2) + \G[\\p^;p,)\\ (42) 

In Eq. (|42p. we have linked G*^^^ with the first order correlation functions of g\^\ The first 
term in Eq. l|42|) is the product of the mean intensities measured by the two detectors. The 
second term, which corresponds to the "intensity fluctuation correlation" in Eq. p5|) . is noth- 
ing but the two-photon interference term. The superposition takes place between quantities 
52(^*21 K)gi{pi, k') and 32 (P2, '5')5i(Pii ^) Eq. (|¥T|) . namely the two-photon amplitudes. 
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Now, we calculate the |G^2''(Pi; P2)P term of Eq. ([55]). Following the experimental setup of 
Fig. [9](b), the Green's function of free-propagation can be written as: 



27rc dB 

where Ps is the transverse vector on the source plane. We have therefore propagated the field 
from the source to the pi plane and p2 plane in arm A and arm B, respectively. 

Substituting k) and 52 (P2, 'S) into (pi; P2) and completing the integral on dn, we 

obtain 

G^i2{Pi'iP2) = / dfi gl{pi,K)g2{p2,K) 

dps e-'^''^ ^-'^\P^-Ps\'' e^^ds gi5^IP2-p-el= (43) 

If we chose the distances from the source to the two detectors equal (d^ — ds), the integral 
of dps in Eq. (j43[) yields a point-point correlation, S{pi — P2), between the pi plane and the 
P2 plane, while taking an infinite size of the source. Thus, after the integration of the bucket 
detector, the joint-detection counting rate is calculated as: 

Rc{P2) cc J dpi \A{p,) S{pi - p2)\^ - \A{p2)\' (44) 

where A{pi) is the aperture function of the object plane. 

Thus we have successfully explained the experimental observation in terms of two-photon 
interference. 

As we know has a counterpart of F^^ in classical optics, namely the first-order coher- 
ence function. Can we use classical interference to explain this phenomenon? In fact, the theory 
of statistical correlation of intensity fluctuations is not a "must" for explaining the second-order 
correlation of thermal radiation. Examining G^^-' — G^'iiG^J + 1^12'' P of thermal light, it is easy 
to find that the intensity fluctuation correlation corresponds to the second term of jG^^ p. G^j"* 
characterizes the first-order coherence, i.e., the ability of observing first-order interference, of 
the fields E{ri,ti) and E{r2,t2). Why don't we explain HBT in terms of classical interference, 
instead of the intensity fluctuation correlation? As we usually say: "people are always smarter 
than theory." This time, however, people have been too smart (we lost 50 years)! There is 
a reason not to do that, because E{ri,ti) and E{r2,t2) would not interfere with each other 
in the case of HBT: £^(ri,ti) and E{r2,t2) are measured by two photodetectors at distance. 
Maxwell's wave equation is a linear differential equation. If E{ri, ti) and E{r2,t2) are solutions 
of the wave equation then the resulting field of E{r,t) — E{ri,ti) + E{r2,t2) is also a solu- 
tion. However, to claim a solution of the wave equation or interference, E{ri,ti) and E{r2,t2) 
have to be physically superposed at a space-time point, as we usually do in an interferometer. 
For example, in the Youngs double slit experiment, E{ri,ti) and E{r2,t2) are the fields at 
the double slits in early times relative to the time of photodetection. The resulting field of 
E{r, t) = E{ri,ti) + E{r2,t2) is physically added at the space-time point of the photodetection. 
In HBT and Scarcelli's experiment, however, the situation is different. E{ri,ti) and E{r2,t2) 
are never brought together; further more, ii^(ri,ii) and E{r2,t2) are measured by two distant 
photodetectors! What kind of interference is this? 

It is interference, but a different type. Under the framework of Glauber's photodetection 
theory, we have provided a two-photon interference picture. It is the interference between two- 
photon amplitudes. Unfortunately, this concept has no counterpart in classical electromagnetic 
theory of light, unless one makes the theory nonlocal by assuming E{ri,ti) and E{r2,t2) are 
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addable nonlocally through the measurement of two photodetectors, or makes a concept such as 
the interference between [field A goes to Di x field B goes to D2] and [field B goes to Di x field 
A goes to D2], which does not make any sense under the framework of Maxwell electromagnetic 
wave theory of light. 

To conclude this section: the two-photon imaging experiment of Scarcelli et al. has demon- 
strated that the ghost imaging of chaotic light is an interference phenomenon involving the 
superposition between indistinguishable two-photon alternatives, rather than statistical cor- 
relation of intensity fluctuations. The two-photon correlations are observed in the intensity 
fluctuations, however, they are not caused by the statistical correlation of the intensity fluctu- 
ations. 



A statement from the author 

This article was originally prepared as lecture notes for my students. The lecture notes were 
turned into a review paper at the end of 2006. The main purpose of this article is to clarify the 
fundamental concepts of quantum imaging and to provide the basic knowledge of physics and 
necessary tools of mathematics, which are involved in quantum imaging, for the general physics 
and engineering community. The debate regarding the quantum or classical nature of quantum 
imaging is still ongoing and will probably endure for a while. In a sense, this is a debate 
that started 50 years ago since the discovery of HBT. In this article, we have emphasized the 
two-photon coherence nature of the quantum imaging phenomena. As usual, there are always 
different opinions J30j. Classical formalisms do exist, such as Eq. ([T]); such formalism, however, 
is far from capturing the essence of the phenomenon and may mislead us from the truth, as 
it has happened to the two-photon physics of thermal light. The interpretation in terms of 
two-photon coherence elegantly explains all of the features of the experiments carried out with 
both "quantum" and "classical" sources. For this reason, two-photon coherence seems more 
suitable for explaining the physics of quantum imaging. 



Appendix A: Fresnel diffraction-propagation 

In Fig. lA— 1[ the field is freely propagated from plane ctq to plane a. For imaging applica- 
tions, it is convenient to describe such a propagation in the form of Eq. ([5]). In this Appendix, 
we evaluate the g{li,uj]p,z), namely the Green's function or the field transfer function, for 
free-space Fresnel propagation. 



E(p,z) 




Figure A—1: Schematic of free-space Fresnel propagation. The complex aperture function 
^(po) is composed by a real function A{pq) and a phase e"*'* ''" associated with each of the 
transverse wavevector k on the plane of Gq. Only one mode of wavevector k(K, oj) is shown 
in the figure. 
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According to the Huygens-Fresnel principle, the field at space-time point (p, z, t) is the 
result of a superposition among these spherical secondary wavelets originated from each point 
on the (To plane, see Fig. \A—1[ 

E[p,z,t) = J dKdujE[K,uo) J e-*('^*-'=''') (A-l) 

where A{po) is the complex aperture function and we have taken zo = and to = as usual. We 
have also assumed a common complex amplitude E(k,, lo) for all point sources of the secondary 
wavelets on the plane of cto- In a paraxial approximation, we take the first-order expansion of 
r' in terms of z and p 

^'^ 
2z2 

E{p, z, t) is thus approximated as 

E{p,z,t)^ f dKdujE{K,to) f dpo^^e'^^^e'^^f^-P^^^e-'"^* {A-2) 



r'^^z^ + \p~po\^:^z{l + \^-^). 



where e^a^'^"^"' is named as the Fresnel phase factor. 

Assuming the complex aperture function is composed by a real function A{po) and a phase 
g-iK-Po associated with the transverse wavevector and the transverse coordinate on the plane 
of ctq, which is reasonable for the setup of Fig. 1[ E{p, z, t) can be written in the following 
form 



E(p,z,t)= J dfidijE(K,uj)e-''^* J dpoA{po)e'^-P° ^-^e'^^P-P''^^ {A-3) 
The Green's function for free-space Fresnel propagation is thus 

gin, u;;p,z)^^— f dpo A{po) e^^^^-P'^^" (A-4) 
where —iLu/2Trc is a normalization constant. 



Appendix B: Quantum state of thermal light 

We assume a large number of atoms that are ready for two-level atomic transition. At most 
times, the atoms are in their ground state. There is, however, a small chance for each atom to 
be excited to a higher energy level and later release a photon during an atomic transition from 
the higher energy level E2 {AE2 7^ 0) back to the ground state Ei. It is reasonable to assume 
that each atomic transition excites the field into the following state: 

|*)^|0)+6^/(k,s)4jO} (B-l) 

where |e| <C 1 is the probability amplitude for the atomic transition. Within the atomic 
transition, /(k, s) = ( \I'k,s | ^I' ) is the probability amplitude for the radiation field to be in the 
single-photon state of wave number k and polarization s: | ^k,s ) — | lk,s ) = oj^. ^ | ). 

For this simplified two-level system, the density matrix that characterizes the state of the 
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radiation field excited by a large number of possible atomic transitions is thus 

p^l[[\0)+eJ2f{k,s) | ) } (B-2) 

toj k,s 

X n { ^ I + ^* E /(k', s') e-'*- ( \au',s' } 

tok k',s' 

^ { 1 ) + E /c^' ^) e^^'^*"^ «L 1 )] + [• ••] } 

X {(0| +e*[^ 5] /(k',s')e^"'*-(Olak',.'] +£*'[...]} 

tofc k',s' 

where 6"*"*"^ is a random phase factor associated with the jth atomic transition. Since |e| ^ 1, 
it is a good approximation to keep the necessary lower-order terms of e in Eq. (|i3— 2D . After 
summing over toj (^ofc) by taking all its possible values, approximately, we have 

p=.|0)(0| + |e|^^|/(k,s)|^|lk,s)(lk,s|, 

k,s 

+ 1^1' E E |/(k',5')|' I lk,slk'.«' )(lk,.lk'..' |. (B-3) 

k.s W .s' 
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